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ABSTRACT 


For  a  communications  system  with  a  low  data  rate  requirement  (i.e.  , 

1  -  1000  bits/sec)  one  can  use  troposcatter  links  over  a  range  interval  beyond 
the  distances  commonly  used.  The  range  is  limited  by  the  requirement  that 
the  tropopause  layer  is  to  be  included  in  the  common  volume  of  the  beams. 
This  effectively  limits  the  range  in  most  cases  to  about  450  miles.  It  is 
proposed  to  use  frequencies  in  the  interval  3  to  5  GHz.  This  enables  one 
to  specify  a  completely  mobile  troposcatter  system  since  the  antenna  para¬ 
boloids  have  diameters  of  the  order  of  10  feet  and  the  power  requirement  is 
a  modest  10  kW. 
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Feasibility  of  Low  Data  Rate  Mobile 
Troposcatter  Communication  Systems 

INTRODUCTION 

During  the  last  decade  the  utilization  of  tropospheric  scatter  communi¬ 
cation  systems  for  military  purposes  has  increased  substantially.  In  general 
these  military  systems  can  be  characterized  by  moderately  large  receiver 
bandwidths  (for  multiple  voice  channels)  and  large  stationary  receiver  and 
transmitter  complexes.  Links  in  the  VHF  band  (50  -100MHz)  provide  tropo- 
scattered  fields  out  typically  to  600  miles  beyond  the  radio-horizon.  Links 
in  the  UHF  and  SHF  bands  (300MHz- 5GHz)  are  typically  used  over  ranges  of 
the  order  of  100  to  300  km.  For  these  large  bandwidth  systems  the  trade-off 
is  primarily  between  larger  antenna  size  and  range. 

However,  if  we  are  allowed  to  restrict  our  receiver  bandwidths  to  a  few 
KHz  we  can  synthesize  a  class  of  troposcatter  communication  systems  which 
have  both  small  antenna  size  (antenna  dish  diameters  of  the  order  of  a  few 
meters)  and  moderately  long  range  capability  (links  of  the  order  of  720  km.  ). 
The  penalty  paid  is  that  one  is  usually  restricted  to  binary  communication 
techniques  with  a  moderately  low  information  bit  rate.  We  propose  to  use  a 
FSK  (frequency  shift  keying)  modulation  scheme.  The  benefits  of  such  a 
troposcatter  link  are  significant.  If  frequencies  in  the  3-6GHz  region  are 
used  the  transmitter  and  receiver  complexes  can  be  made  mobile  because 
of  the  small  antennas  and  modest  transmitter  power  requirements.  Moreover 
at  distances  of  the  order  of  640  km,  the  lower  boundary  of  the  common 
volume  of  the  intersecting  transmitter  and  receiver  antenna  beams  is 
typically  several  kilometers  above  the  earth's  surface.  The  scattering 
region  is  then  in  the  vicinity  of  5  to  20  kilometers  of  altitude.  Accordingly, 
the  communications  link  will  not  be  degraded  by  an  adverse  condition  of  the 
ionosphere.  Finally,  hardware  needed  for  such  a  system  is  both  relatively 
inexpensive  and  well  within  the  present  technological  capability. 
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This  report  is  divided  into  two  parts.  Part  I  is  in  part  tutorial  describ¬ 
ing  the  nature  of  the  atmosphere  which  is  pertinent  to  our  proposed  tropo- 
scatter  link.  The  effects  of  refraction  are  accounted  for  and  diffraction  is 
shown  to  be  negligible.  For  computational  purposes  a  hypothetical  tropo- 
scatter  link  with  a  communications  distance  requirement  of  400  miles  (640 
km)  is  selected. 

Part  II  is  concerned  with  the  computation  of  the  required  transmitter 
power  levels.  In  the  various  sections,  sufficiently  detailed  analyses  of  the 
contributing  terms  are  given  permitting  an  eventual  parameter  study  to 
decide  on  an  'optimum1  system.  The  bulk  of  the  report  is  concerned  with 
this  area. 
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I.  Qualitative  Look  at  the  Actual  Troposcatter  Geometry  and  the  Common 

Volume  Environment 

In  troposcatter  systems  analysis,  the  appropriate  propagation 
equations  (see  Appendix  A)  are  dependent  only  on  the  fine  scale  variations  of 
the  permittivity  (or  refractive  index).  However,  these  fine  scale  fluctuations 
are  dependent  on  the  gross  meteorological  conditions.  The  large  scale 
refractive  index  behavior  is  described  in  Section  A  for  its  tutorial  value. 

Then  an  approximation  is  used  for  dN/dh  in  order  to  compute  the  effects  of 
normal  refraction.  By  a  suitable  transformation  the  rays  are  "straightened  ,r 
in  Section  B.  Finally  the  beam  geometry  in  relation  to  the  actual  path  can  be 
presented.  A  method  for  computing  the  scattering  angle  ft  is  given. 

If  the  reader  is  interested  only  in  the  system  sizing  considerations,  he 
is  advised  to  proceed  directly  to  part  II  of  this  report. 

A.  Atmospheric  Structure 

In  the  design  of  transhorizon  communication  links,  a  relatively 
detailed  knowledge  of  the  atmospheric  structure  is  essential.  In  particular 
the  volume  of  atmosphere  delineated  by  the  intersection  of  the  transmitter 
and  receiver  main  lobe  antenna  beams  is  of  paramount  importance.  In  order 
to  quantitatively  describe  the  energy  "scattered"  out  of  this  common  volume 
it  is  necessary  to  estimate  the  structure  of  the  refractive  index  variations 
both  in  time  and  space.  This  is  necessary  in  order  to  catagorize  whether 
the  "scattering  "  mechanism  is  representative  of  specular  reflection,  diffuse 
reflection  or  the  usual  inhomogeneity  scattering  (or  some  combination  of 
these  processes). 

It  will  be  seen  that  for  a  400  mile  link  the  common  volume  will  straddle 
the  upper  troposphere  and  the  lower  stratosphere  (see  Fig.  1).  It  is  expected 
that  refractive  index  fluctuations  in  the  tropopause  vicinity  will  give  a  signifi¬ 
cant  received  signal.  This  is  in  contrast  to  the  usual  shorter  transhorizon 
links  where  most  of  the  energy  is  scattered  at  altitudes  below  5  km.  Because 
the  common  volume  of  a  400  mile  link  is  so  high,  very  little  of  the  vast 
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Fig.  1.  The  U.S.  extension  to  the  I.C.A.O.  standard  atmosphere. 
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amounts  of  experimental  troposcatter  data  collected  over  the  past  decade  is 
directly  applicable.  However,  we  can  still  size  the  system  based  on  the 
limited  experimental  data  available.  If  better  estimates  are  desired  it 
would  be  justified  to  perform  some  basic  experimental  research  for  links  of 
this  distance  and  operational  frequency. 

The  refractive  index  structure  will  specify  the  type  of  radio  wave 
propagation.  For  practical  radio  scatter  applications  the  refractive  index 
is  usually  described  in  terms  of  N  units  (rather  than  the  refractive  index  p  ) 
where  N  is  defined  by  [  1]  , 

N  =  (v  -  1)10+6  =-  77.6  ^  -  5.6-y  +  3.75  x  105-^-  (1) 

where  e  is  the  partial  pressure  of  water  vapor  (mbars) 

P  is  the  total  pressure 
T  is  the  temperature. 

To  a  first  approximation 

N  -  77.6  y  +  3.75  x  105  =  D  +  W  (2) 

where  D  refers  to  the  dry  term  and  W  to  the  wet  term.  Figure  2  [  1]  shows 
that  only  the  dry  term  is  important  above  an  altitude  of  5  km.  Of  course, 
it  is  only  the  fluctuations  in  N  (or  D)  that  is  of  interest.  Unfortunately,  not 
many  refractivity  profiles  in  the  upper  troposphere  and  lower  stratosphere 
are  available  by  direct  measurements  (particularly  the  fine  scale  fluctuations). 
This  being  the  case  it  is  best  to  choose  turbulence  models  which  seem 
applicable  at  lower  altitudes  and  check  whether  the  limited  long  link  experi¬ 
mental  scatter  data  is  consistent.  It  must  be  realized  that  the  performance 
of  a  communications  link  is  dependent  on  the  gross  (as  compared  to  the  fine 
scale  turbulence)  meteorological  conditions  which  vary  both  diurnally  and 
seasonally.  For  example  there  is  evidence  that  strong  scattering  occurs  in 

'  In  short  link  troposcatter  links  it  is  the  fluctuations  in  W  that  are  dominant. 
For  the  longest  links  the  D  term  is  dominant  and  should  be  relatively  insensi¬ 
tive  tc  the  actual  path  chosen. 
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the  tropcpause  region.  This  layer  shows  a  marked  drop  in  altitude  during 
the  winter  and  a  rise  in  the  summer.  This  points  out  although  we  will  use 
a  rather  simple  turbulence  model  in  this  report,  any  more  extensive  link 
analysis  must  consider  the  larger  scale  meteorology  (in  particular  the  stable 
and  unstable  layer  distribution). 

B.  Refraction  and  Diffraction 

Since  the  atmospheric  density  decreases  systematically  with 
altitude,  an  oblique  radio  wave  will  undergo  refraction.  The  local  curvature 
of  the  radio  ray  is  given  by 

c  =  cos  e  (3) 

where  7;  and  ^  are  local  values  of  the  refractive  index  and  the  elevation 
angle.  To  simplify  the  geometry  of  the  propagation  problem  it  is  customary 
to  straighten  out  the  actual  curved  path  of  the  radio  ray  by  presenting  the 
straight  ray  relative  to  an  imaginary  earth  of  radius  k  r^  where  the  relative 
curvature  between  the  ray  and  earth  is  maintained 


r 

e 


dr? 

dh 


cos 


1 

kr 

e 


+  0 


earth  ray  effective  earth 

curvature  curvature  curvature 


curvature  of 
straight  ray 


If  the  "standard  atmosphere"  gradient 


(4) 


d72_  _  _  1 

dh  4  r 

e 

is  used  then  k  =  4/3  and  the  effective  earth  radius  is  4/3  r^  =-  8500  km. 
The  atmosphere  in  this  imaginary  earth  will  have  a  constant  mean  77  with 
altitude.  Figure  3  [2]  shows  that  for  the  altitude  regime  0  to  17  km,  the 
(4/3  r^)  approximation  gives  values  of  refraction  which  are  less  than  the 
actual  average  bending.  Retaining  the  (4/3  r  )  approximation  will  then  give 
us  a  conservative  estimate  which  is  a  desirable  feature. 
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Fig.  3.  Refraction  cor rections  . 
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Two  types  of  diffraction  occur  in  transhorizon  propagation;  namely 
diffraction  around  a  smooth  spherical  surface  and  knife  edge  diffraction. 
For  a  smooth  surface  the  distance  beyond  the  radio  horizon  at  which  the 
scattered  energy  and  diffracted  energy  are  approximately  equal  is  given  by 


d 


km  where  f  is  in  GHz. 


2  1  km. 


For  f  -  3  GHz 

(5) 


Since  the  horizon  distance  is  of  the  order  of  10  km,  it  can  be  safely  assumed 
that  the  scattered  fields  will  dominate  at  distances  greater  than  31  km  from 
the  transmitter.  At  a  few  hundred  kilometers  the  diffracted  fields  are 
entirely  negligible.  Knife  edge  diffraction  would  be  significant  only  if  the 
lower  edge  of  the  common  volume  is  directly  at  or  below  the  top  of  some 
obstacle.  The  section  on  beam  geometry  will  show  that  this  is  not  a  detect¬ 
able  propagation  component. 

C.  Path  Considerations 


It  is  a  well  known  experimental  fact  that  the  received  scattered 
power  in  a  troposcatter  link  decreases  rapidly  with  an  increasing  scattering 
angle  g  (ft  is  defined  on  Fig.  4).  The  reason  for  this  is  that  not  only  are 
the  theoretical  scattering  models  very  ft  dependent  but  also  the  intensity  of 
the  turbulence  tends  to  decrease  in  magnitude  for  increasingly  higher  common 
volumes  (that  is  AN^  is  altitude  dependent  which  for  a  bistatic  link  is 
impossible  to  separate  from  the  ft  dependence).  The  objective  then  is  to 
keep  the  beam  as  close  to  the  local  horizon  as  possible  and  at  the  same  time 
guard  against  obstacle  blockage  of  the  antenna  beam. 

To  see  the  effect  that  link  distance  and  obstacles  have  on  ft  let  us 
first  consider  a  simplified  geometry.  Assume  that  the  path  has  a  uniform 
surface  elevation  of  500  ft.  above  sea  level  and  moreover  that  the  antenna  is 
placed  2D  ft.  above  this  surface  (different  antenna  heights  or  surface  levels 
can  be  easily  accounted  for  in  most  cases  by  a  radial  translation).  Define 
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Fig.  4.  Link  geometry. 
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the  horizon  line  to  be  that  line  passing  through  the  antenna  dish  center  and 
being  tangent  to  the  surface.  The  elevation  angles  ^  are  measured  with 
respect  to  this  line.  Figure  5  shows  the  maximum  permissible  height  of 
obstacles  as  a  function  both  of  elevation  angle  ^  and  distance  from  the  trans¬ 
mitter  (cr  receiver).  Figure  6  shows  a  similar  set  of  curves  which  gives 
the  height  of  the  common  volume  provided  both  the  common  volume-receiver 
(or  transmitter)  distance  and  elevation  angle  ^  is  specified.  Here  ^  is 
understood  to  be  the  elevation  of  the  bottom  edge  of  the  main  antenna  lobe. 
The  scattering  angle  g  can  be  computed  from  the  following  set  of  equations 


8  -  m )  f 

”  re 


(6) 


.  -1  ,re  +  1t  ,  .  .  -l,re+1r  , 

5t =  sin  ■  6r  =  <r-nr> 


(7) 


6  n 

=  fl  +  (6t+et+  j-1  -  )  +  (6  r  +  e  r  +  -J’  r  )  -  180 


(8) 


For  the  case  where  the  horizon  is  several  km  away  6  ^  =“  6  =**  90°  and 


0V  t  r 

B  ~  d  +  (€t  +  -f-)  +  (cr  +  -f~) 


(9) 


The  normal  objective  is  to  make  (?r+  as  small  as  possible  so  that  the 

tropopause  is  included  in  the  common  volume.  Communication  across 

mountain  ranges  is  easily  possible  provided  the  receiver  and  transmitter 

*  locations  are  positioned  sufficiently  far  enough  away.  The  standoff  distances 

can  be  determined  from  Figs.  5  and  6.  For  most  areas  of  the  United  States 

it  is  easy  to  find  terminal  locations  where  both  e  and  e.  can  be  made  1/4° 

r  t 

or  less.  Accordingly  cr  =  =  l/4°  will  be  used  in  the  calculations.  For 

a  tropo scatter  link  of  640  km,  the  lower  edge  of  the  common  volume  will 
then  be  at  an  approximate  altitude  of  7  km.  This  is  sufficiently  low  such 
that  the  tropopause  layer  will  be  included  in  the  common  volume  for  the 
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-PERMISSIBLE  OBSTACLE  HEIGHT  ABOVE  SEA  LEVEL  (km) 


d0  -  DISTANCE  FROM  OBSTACLE  TO  RECEIVER  (or  Transmitter)  ( km) 


Fig.  5.  Permissible  obstacle  height  versus  distance  from 
receiver  (or  transmitter). 
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FEET 


Fig.  6.  Common  volume  altitude  . 
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entire  year.  The  parameters  to  be  used  in  the  sample  calculations  are  then 


d 

640  km 


X 

320  (km) 


ft  =  cr 


0.25 


e  =e 

V,  t  V, 

1.5° 


/3 

~6.  2° 
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II.  System  Sizing 


A.  Method  of  Calculation  of  the  Required  Transmitter  Power 

One  purpose  of  the  report  is  to  provide  enough  information  on  the 
hardware  and  propagation  parameters  of  a  (400  mile)  troposcatter  communi¬ 
cations  link  so  that  one  can  make  quick  computations  for  alternative  systems. 
The  selection  of  a  system  may  be  constrained  by  a  number  of  factors  external 
to  what  we  have  considered  and  consequently  we  have  not  done  an  all  encom¬ 
passing  parameter  study.  What  we  intend  to  do  is  to  choose  an  acceptable 
receiver  system,  antenna  size  and  beam  geometry  and  to  evaluate  the  average 
power  required  for  the  transmitters.  If  these  power  levels  are  sufficiently 
low  (say  10KW  or  less)  the  system  will  be  considered  feasible.  In  computing 
the  required  transmitter  power  level  we  will  account  for  the  various  quantities 
listed  in  table  1  (and  discussed  in  detail  in  the  identified  section). 

The  formula  for  minimum  transmitter  power  is  listed  below.  The 
identification  and  magnitude  of  the  terms  in  the  formula  are  discussed  in  the 
appropriate  sections  (all  terms  in  dB).  The  signs  associated  with  the  terms 
below  are  consistent  with  the  signs  on  the  figures  in  the  appropriate  sections 


P 


T 


Minimum  Effective  Antenna  Gain 

Received 

Signal 


Free  Rain  Bistatic 

Space  &  Cross 

Loss  Gas  Section 


Atten. 


(10) 


Distance 
Factor  of 
Bistatic 
Radar 


A  typical  computation  will  be  shown  in  section  II- 1. 
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TABLE  I 

Quantity  Section 

Receiver  and  antenna  noise  -N  =  kTQp  II- G 

Antenna  gains  II- B 

Antenna  efficiency  L^  =  Lp  =  3dB 

Aperature  -  to  -  medium  coupling  L  II-C 

a 

Rein  and  Gaseous  Attenuation  L  -  II-D 

K ,  O 

Free  Space  Attenuation  Lp^  II-E 

Scattered  Power  normalized  P/Pfs  II-F 

by  the  Free  Space  Power 

Signal  to  Noise  Ratio  P  /N  II-H 

nr  o 
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B. 


Antenna  Considerations 


Because  of  the  requirements  of  mobility  it  is  necessary  that  the 
antennas  be  relatively  small.  If  the  antenna  is  say  less  than  5  meters  in 
diameter  one  can  easily  envision  systems  where  any  or  all  of  the  transmitters 
and  receivers  can  be  made  mobile  by  mounting  on  trucks.  Since  troposcatter 
links  require  beamwidths  of  the  order  of  a  degree,  this  requires  that  the 
frequency  of  operation  be  in  the  GHz  region.  Paraboloids  are  efficient  reflectors 
at  these  frequencies  and  are  the  only  antennas  considered.  The  frequency 
range  of  interest  is  from  1  GHz  to  8  GHz  (the  lower  limit  is  constrained  by 
the  antenna  size  and  the  upper  limit  by  the  attenuation  due  to  rain). 

In  free  space  the  appropriate  gain  expressions  are  given  by 

(a)  the  directive  gain  G^ 


G 


max  rad.  intensity 
D  ~  average  rack  intensity 


If  0  ,  0  are  the  half-power  beam- 

H  V 


widths  in  two  orthogonal  plans  the  directive  gain  can  be  approximated  by 


~  41253  /0  q  j  x 

GD=  TOTT  (0H’  9V  are  in  desrees) 
H  v 


on 


(b)  the  power  gain  G 

_  max  rad.  intensity 

rad.  int.  from  lossless  isotropic  Source  with  same  power  input 

(c)  The  radiation  efficiency  factor  p 


<  1 


(12) 


(d) 


Effective  aperature  and  antenna  aperature  efficiency 
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A  is  physical  size  of  antenna.  (13) 


4rr  A 


4np  A 


Combining  equations  11,  1Z  and  13  gives 


Vv  =  41253 


XZD 

r 

4tt  Ao” 


For  the  paraboloid  of  diameter  D 


the  approximate  mainlobe  beamwidth  is 

9  =-7^7  =  ^53^(]-d)  [degrees]  (14) 

Figure  7  gives  the  beamwidth  as  a  function  of  frequency  with  the  diameter 
of  the  paraboloid  (D)  as  a  parameter.  The  antenna  gain  divided  by  the 
aperature  efficiency  is  shown  in  Fig.  8. 

C.  Aperture-  to-Medium  Coupling  Loss 

The  theoretical  antenna  gain  expressions  just  described  presuppose 
that  the  incident  wave  is  an  ordinary  plane  wave.  However,  in  scatter  propa¬ 
gation  the  individuals  scatters  in  the  common  volume  produce  a  wave  at  the 
receiving  antenna  which  is  not  planar  and  hence  the  full  gain  of  the  antenna 
will  not  be  realized.  The  incoherence  of  the  wave  front  causes  an  effective 
gain  reduction  which  is  greater  for  the  large  antenna  than  for  the  small 
antenna.  There  exists  no  simple  relationship  which  accurately  gives  the 
coupling  loss  in  terms  of  the  antenna  parameters  and  geometry  alone.  Hence 
we  must  estimate  this  loss  by  comparison  of  the  proposed  antenna  parameters 
with  existing  experimental  estimates  of  aperture -to-medium  coupling  losses. 
Yeh  [3]  has  collected  a  large  number  of  experimental  estimates  and  plotted 

-zr-  (Fig.  9)  For  a  typical  case  in  the 

GHz  region  for  the  400  mile  link  we  have 


them  versus  the  parameter  v>/ 
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MAIN  LOBE  BEAMWIDTH  (deg) 


Fig.  7.  Beamwidth  versus  frequency. 


19 


Fig.  8.  Antenna  gain  versus  frequency. 
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Lo(dB)  —  APERTURE -TO -MEDIUM  COUPLING  LOSS 


Fig.  9.  Aperture -to -medium  coupling  loss  . 
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i _ =  z° _  >  4.1 

jrr.  5)u(i.  5u) 


This  gives  an  aperture -to- medium  coupling  loss  of  about  9.  OdB  which  is 
probably  slightly  on  the  conservative  side. 

D.  Rain  and  Gaseous  Absorption  Loss 


In  the  GHz  region, attenuation  due  to  rain  and  the  gaseous  absorp¬ 
tion  can  become  significant.  Fig.  10(Bean  [  1]  )  shows  the  total  attenuation 
due  to  these  two  sources  as  a  function  of  frequency.  The  magnitude  of  the 
values  of  attenuation  will  be  less  than  the  curves  99$  of  the  time.  The 
400  mile  values  can  be  easily  interpolated  from  the  two  curves  for  300  mile 
and  1000  mile  path  lengths.  For  example  at  4GHz  the  total  path  attenuation 
is  about  4.  5  dB.  The  rapidly  increasing  attenuation  with  frequency  is  the 
fact  that  restricted  our  investigation  to  frequencies  less  than  8GHz. 

E.  Free-Space  Attenuation 

In  scatter  theory  the  received  power  is  usually  normalized  by 
the  power  that  would  be  received  over  a  line-of- sight  path.  For  lossless 
propagation  the  free  space  power  is  given  by 


Pfs  PTGTGR 


1 6n^d 


(15) 


The  basic  propagation  attenuation  in  free  space  is  defined  by 

P, 


'Pf 


f  s 

ptgtgr 


-  (Jw)2 

4nd 


(16) 


L  ^  is  plotted  in  Fig.  11  with  d  as  a  parameter.  Now  the  theoretical  gain 
product  should  be  modified  by  the  aperture -to-medium  coupling  loss  La. 
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FREQUENCY  (MHz) 


Fig.  10.  Combined  rain  and  gaseous  absorption  (l%exceedance  values). 
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Fig.  11.  Free  space  attenuation. 
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The  free  space  power  is  then  given  by 


Pfs[dB]  =  PT  -  Lpf  +  |GT+GR-La|  (all  terms  in  dB)  (17) 

F.  Scattering  Theory 

In  a  previous  section  we  have  accounted  for  refractive  corrections 
and  excluded  diffraction  as  being  important. 'Scattering1  is  often  used  as  an 
alt  inclusive  name  in  transhorizon  propagation  to  describe  any  energy  which 
leaves  the  common  volume  region  in  a  direction  different  from  the  axis  of 
the  transmitter  beam.  For  long  distance  links  we  are  fortunate  in  that  only 
one  mode  (the  usual  inhomogeity  scattering)dominate s. 

Various  physical  models  for  transhorizon  propagation  were  developed 
during  the  1950's.  [5]  During  that  time  there  was  major  debate  on  whether 
the  appropriate  mechanism  was  specular  reflection  from  elevated  layers  or 
scattering  from  the  refractive  index  fluctuations  caused  by  the  atmospheric 
turbulance.  Gradually  the  adherents  of  the  scattering  theory  became  more 
prevalent  even  though  no  single  mathematical  scattering  model  consistently 
agreed  with  the  experimental  results.  In  the  early  sixties  extensive  experi¬ 
mental  work  by  the  French  (see  DuCastel  [4])  revealed  that  the  received 
signal  contained  a  mixture  of  scattered  in  addition  to  diffuse  and  specularly 
reflected  energy.  Since  then  the  main  effort  has  been  to  obtain  a  viable 
syntheses  of  these  diverse  theories.  For  long  distance  propagation  (if  we 
neglect  anomalous  effects  which  are  not  too  uncommon)  the  recent  tendency 
has  been  to  describe  the  physical  process  in  terms  of  a  standard  scattering 
model  and  to  include  the  effects  of  any  layering  in  the  autocorrelation  coeffi¬ 
cient  of  the  refractive  index  variations.  This  will  be  the  approach  used  here. 

The  fact  that  energy  is  reflected/scattered  from  layered  regions  is  now 
well  established.  For  long  distance  transhorizon  propagation  the  layers  at 
the  tropopause  level  provide  a  good  scattering  region  (this  is  certainly  true 
for  the  400  mile  links  considered  in  this  report.  )  To  a  lesser  extent  layers 
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in  the  stratosphere  also  contribute  to  the  scattered  return. 

In  Appendix  A  we  have  outlined  the  derivation  of  the  equation  for  the 
average  received  power  normalized  by  the  free  space  received  power  (equation 
a-20).  The  electrodynamics  in  this  derivation  are  complex  but  are  relatively 
straight  forward  requiring  few  physical  assumptions.  We  have,  however, 
been  forced  to  assume  local  homogeneity  in  the  statistics  of  the  refractive 
index  fluctuations.  This  barely  tolerable  assumption  is  made  because  the 
alternative  leads  to  a  prohibitively  difficult  problem.  If  the  beams  are  very 
narrow  this  is  however  a  reasonable  assumption  to  make.  The  remaining 
difficulty  is  to  specify  the  form  of  the  spectrum  (<f>n  (K))  of  the  refractive 
index  fluctuations.  Many  investigators  over  the  years  have  postulated 
various  mathematical  models  for  (K).  Usually  these  models  took  the 
form  of  (K)  ~  K  m  where  m  is  a  fixed  constant.  However  none  of  these 
models  have  been  universally  accepted  because  the  experimental  evidence 
often  seemed  contradictory.  It  was  soon  realized  that  for  physical  reasons 
the  refractive  index  spectrum  had  a  variable  wavelength  dependence,  f  6] 

Fig.  12  shows  such  a  spectrum  broken  up  into  three  regions  each  with  a 
different  exponential  dependence.  At  present  the  most  widely  used  exponent 
is  m  =  -**/3  which  refers  to  the  inertial  subrange.  It  should  be  noted  that  this 
dependence  was  first  arrived  at  by  Kolmogorov  in  1941  on  the  basis  of  a  simple 
similarity  analysis  of  the  turbulence,  (see  Tartarski  [7]  and  Bachelor  [8] 
for  details).  In  most  of  the  cases  where  experimental  evidence  seems  to 
show  that  the  ,T  1  1  / 3 ”  dependence  law  "fails"  it  is  often  because  the  value  of 

K(=  2k  sin  3/2)  for  that  particular  experiment  lies  outside  the  inertial  subrange. 

*  2 

Inasmuch  as  we  have  assumed  that  is  dependent  on  position  we  will 

plot  f  <fr(K,  y)dy  rather  than  <f>(K,  y)  since  there  exists  numerical  values  derived 
Vy 

for  the  former  but  not  the  latter  quantity. 
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Fig.  12.  Log1()  /  4>  (K,  7)  dy  versus  scale  size  I . 
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This  is  particularly  true  for  the  longer  wavelength  and  shorter  path  links 
where  both  k  and  0  are  smaller  (this  means  that  the  characteristic  eddy 

size  1  =  is  typical  of  the  buoyancy  subrange  which  has  a  different  value 

of  n).  Most  of  the  common  troposcatter  links  fall  in  this  category  and  hence 
this  is  to  a  large  part  the  reason  why  the  standard  scattering  theories  were 
not  well  accepted.  For  the  GHz  range  we  are  usually  in  the  inertial  range. 
For  example  if  we  take  f  =  3GHz  and  0  =  6.2°  the  characteristic  (scale) 
size  is 


1  =£  •  2  sin  3/2  -  -^meters  (]8) 

This  value  of  1  is  well  within  the  inertial  subrange  and  hence  the  11 1  1  /3 11  law 
should  b£  appropriate.  Tartarski  [7]  gives  the  spectrum  as 

<f>  (K)  =  (.  033)  C2  K'11/3  (19) 

n  n 

When  this  is  substituted  into  equation  a-21  the  bistatic  scattering  cross  section 
is 

a  =  (.  384)(sin2X)\"  1/3(sin|)‘11//3  £  C2(y)d  y  (20) 

Y 

The  geometry  and  frequency  being  specified  it  only  remains  to  determine  the 

volume  integral  f  C  (y)dy  .  This  of  course  can  only  be  determined  experi- 
Vy  n  —  -> 

mentally. 

Crane  [9]  has  recently  undertaken  an  experiment  which  has  given  some 

quantitative  estimates  to  the  volume  integral  C  (^)d^  •  The  entire  experiment 

Vy  n 

will  be  first  described  in  moderate  detail  since  the  parameters  are  very  similar 
to  the  scatter  links  we  are  proposing  and  accordingly  provides  the  feasibility 
for  undertaking  our  troposcatter  link.  His  experiment  utilizes  both  monostatic 
and  bistatic  scattering  techniques.  The  narrow  beam  Millstone  L  band  radar 
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was  used  to  determine  the  spatial  distribution  of  the  layers  through  back 
scattering  cross  section  measurements.  Then  an  X-band  forward  scatter 
link  (between  Wallops  Island,  Va.  and  Westford,  Mass.  )  provided  the 
normal  transhorizon  data.  The  important  parameters  of  the  X-band  forward 
scatter  link  are 

(a)  frequency  ^  7.  7GHz 

(b)  c.  w.  power  =*  lkw 

(c)  minimum  detectable  signal  =  -190dBw  (signal  levels  ranged  from 
-174dBw  to  -190dBw) 

(d)  transmitter  antenna  has  6ft.  diameter  and  varied  in  elevation 
while  the  receiver  antenna  had  a  60ft.  diameter  with  a  (.  8°)  elevation  angle. 

(e)  path  length  =  628km  (390  miles) 

The  significant  results  of  the  experiment  were 

(a)  Layers  were  always  detected  near  (and  above)  the  tropopause 

(b)  This  layer  showed  a  relatively  uniform  backscattering  cross 
section 

(c)  Turbulent  layers  were  detected  up  to  22km  (well  into  the  strato¬ 
sphere) 

(d)  The  signal  in  the  forward  scatter  link  in  terms  of  its  scattering 
angle  dependence  was  consistent  with  the  well  known  "  1  I / 3 M  law  of  turbulent 
scattering. 

2 

(e)  From  radiosonde  data  and  the  radar  data,  estimates  of  C  for 

v  '  n 

the  layers  were  made.  Because  of  the  resolution  of  the  L  band  pulse  length 
only  (C^  1^)  could  be  determined  (where  1  is  the  layer  thickness).  A  typical 
value  of  (C^  1^)  is  10  ^(m*^).  Then  ^  (y)dy^  1  0  *  [  m]  ^ ^  where  A^ 

Y 

is  the  horizontal  layer  area  limited  by  the  common  volume. 

If  we  consider  the  geometry  appropiate  to  our  problem,  the  area  limited 
by  the  common  volume  (assuming  that  the  tropopause  is  at  10km)  is  given  by 
(see  Fig.  19  in  Appendix  B) 
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A  _  (AZ>D^h  ^  (10-6.  5)650(1.  5)n  ,  n6 

A  '  Z~n(J/l)  2  tan  (Z.  7) - (7 W  X 

and  hence  the  volume  integral  is  computed  as 

^  Cn  (y)  ~  6.  3  X  10  4  [m7/3] 

Y 

Therefore  from  equation  20  we  have 


6.  26 x  lO8 


2 

m 


a  =*  (2.41  x  10  4)  (sin  \)\  1///3(sin-| )”  1  (21) 

We  will  ignore  the  small  polarization  effect  (since  82°<  y^^90°)  and  compute 
a  as  a  function  of  \  for  a  few  different  values  of  p.  This  is  seen  in  Fig.  13 

We  are  attempting  to  compute  P(d)/P^g.  This  can  be  computed  with 
the  use  of  equation  a-22. 


™/Pfs  =  a[  dB]  +  10  log  - - - 7 - 7 

fS  10  (4tt)  lx"  |Z|d-X'  |Z 

If  we  assume  d=650km  and  lx  1 1  =  ld-xr  |  =  325km  then 

,2 


(22) 


1  0  log 


i  |2 


(4rr)  \X*  r|d-x'  I 


=  -  1  15.  2dB 


and 


P(d)/Pf  =  a  -1  15.  2 

As  an  example  at  f  =  3GHz,  |3  =  6.  2°  Fig.  13  gives  a  value  of  13.  6dB  and 
hence  P(d)/P^s  =  -101.  6dB.  Figure  14  shows  some  experimental  results 

collected  for  different  links  by  N.  B.  S.  Our  value  of  P(d)/P^g  would  seem 
to  be  consistent  with  this  data. 
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BISTATIC  CROSS  SECTION  a  (dB) 


Fig.  13.  Scattering  cross  section  a  versus  frequency. 
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Fig.  14.  Spread  of  experimental  results. 
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G.  Noise  Considerations 

Noise  sources  are  often  classified  as  originating  extra-terrestially, 
terrestially  and  internally  to  the  receiver  system.  In  the  GHz  region  cosmic 
noise  is  insignificant  as  compared  to  the  receiver  noise.  Likewise  man-made 
and  lightning  noise  can  be  ignored  in  this  frequency  regime.  Unfortunately 
we  cannot  ignore  the  earth’s  radiation  since  the  main  lobe  is  almost  tangential 
to  the  horizon  and  the  side  lobes  illuminate  the  earth’s  surface.  Moreover 
because  of  the  effects  of  aperature  blockage  by  the  antenna  feeds  it  is  seldom 
possible  to  get  the  nearest  sidelobes  down  more  than  15dB  relative  to  the 
main  lobe.  The  side  lobes  illuminating  the  earth  (whose  —  300°K)  then 
contribute  to  the  effective  antenna  temperature. 

Let  us  estimate  the  effective  antenna  temperature  for  the  proposed 
beam  geometry  (e~0°.  The  effective  antenna  temperature  is  computed 
from 


T 

a 


JJ  G(0,cd)T^(9,  cd)  sin  0  d0  dcp 


(23) 


where  T^  is  the  brightness  temperature  of  the  Rayleigh- Jeans  law.  From 
Fig.  15  it  is  seen  that  when  the  antenna  beam  is  tangent  to  the  earth’s 
surface  the  brightness  temperature  in  this  direction  is  about  120°K  (for 
the  GHz  region  considered.  Now  if  we  are  using  narrow  beam  high  gain 
antennas  (8~1.  5°,  Gmax^lO,  000),  then  the  average  gain  in  directions 
other  than  the  main  lobe  can  be  crudely  approximated  by 


G 


AV~ 


2 

- zr 

1+cos  /2 


G  max(  1  -  cos  ^/Z) 
(1  +COS  /Z) 


(■  56) 


Since  about  1  /Z  of  the  sidelobe  structure  looks  at  the  ground  (at  T^.  =  300°K) 

and  the  other  l/Z  looks  at  the  sky  (TCI.V  ...  ~  100°K)  we  can  get  an  approximate 

dL.  Y ,  A  V 

antenna  temperature  of 


33 


EQUIVALENT  NOISE  TEMPERATURE  (°K) 


Fig.  15.  Noise  levels. 
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T 

a 


1 


gavtc 


+  2tt 


(,AV  TSKY,  AV 


+  Gmax  (rr  ('5/2)Z)THOR  ^ 


(.  28)  T„  +  (.  28)  T 


SKY,  AV 


+  .  44  T 


HOR 


(24) 


-  (.  28)  (300°)  +  (.  28)  (100°)  +  (.  44)  (120°) 

=-  164°K 

Since  the  antenna  temperature  is  relatively  high  (by  radio  astronomy 
standards)  it  is  not  worthwhile  employing  an  ultra  low  noise  receiver  system. 
Instead  we  choose  as  the  receiver  an  uncooled  parametric  amplifier  with 
an  effective  noise  temperature  of  200°K.  The  total  operating  noise  temperature 
i 5  then 


T 

op 


+  T 


eff 


364° 


The  noise  power  per  1  Hz  bandwidth  is  then  given  by 


N  =  k  T  and  in  terms  of  dB  this  is  (25) 

op 

o  /  4O 

N  =-204  +  1 0  log  n(— — —  )  =“  -203  dB/Hz 
1U  290° 

The  next  factor  one  must  consider  is  the  bandwidth  of  the  individual  receiver 
channels.  The  required  bandwidth  is  dependent  on  several  factors  which  will 
be  examined  in  the  following  section. 

H.  Signal  to  Noise  Ratio  and  Bandwidth  Considerations 

There  are  many  factors  which  must  be  taken  into  account  before 
specifying  the  signal  to  noise  ratio  which  must  be  known  in  order  to  determine 
the  transmitter  power  level.  Let  us  investigate  in  turn  the  transmitter,  the 
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propagation  medium  and  the  receiver. 


Transmitter  -  It  is  initially  assumed  that  the  bandwidth  for  each  of  the 

receiver  channels  is  3KC  (this  will  be  justified  later  in  the  section).  The 

transmitter  will  operate  in  a  FSK  mode  (with  2  n^  possible  frequencies, 

n^  being  the  number  of  chips  per  bit  of  information).  Transmitter  tubes  which 

are  self  oscillators  (like  the  magnetron)  usually  do  not  have  the  long  term 

or  short  term  frequency  stability  necessary  to  keep  the  transmitted  frequency 

within  the  required  passband  of  the  individual  receivers.  Accordingly  we 

specify  that  a  klystron  with  external  (crystal  controlled)  oscillator  with 

stabilities  of  1  part  in  10  be  used.  This  gives  a  maximum  frequency  deviation 

~  ±30Hz  at  3GHz.  (2n  such  oscillators  will  be  needed  in  the  transmitter  as 

P 

well  as  in  the  receiver). 


Propagation  Path  -  For  the  path  geometry  considered  the  following  facts 
are  cogent 

(a)  For  path  distances  greater  than  300  miles,  the  short  term  statistics 
(less  than  1  hr.  )  reveal  that  the  received  signal  envelope  is  very  close  to  being 
Rayleigh  distributed  and  the  phase  has  a  uniform  distribution. 

(b)  The  differences  in  path  length  for  any  two  rays  is  extremely 
small.  For  example  with  6.  2°,  0  ~  1.5°  and  D  =  650  km  equation 

b-4  gives  a  path  difference  of  .  92  km  which  corresponds  to  a  3|jsec  delay. 

This  means  that  the  receiver  integration  times  can  safely  be  brought  down 
to  less  than  a  tenth  of  a  millisecond  if  desired.  The  channel  bandwidth  depends 
indirectly  on  these  differences  in  path  length.  The  small  delay  suggests  the 
bandwidth  may  be  of  the  o'rder  of  a  MHz.  (Unfortunately  experimental  attempts 
to  determine  the  bandwidth  of  tropospheric  channels  have  been  very  sparse  and 
experimentation  in  this  area  should  be  undertaken.  ) 


(c)  The  doppler  frequency  due  to  motion  of  the  atmospheric  irregularities 
cause  band  spreading  at  most  of  only  of  the  order  of  tens  of  Hz  and  accordingly 
can  be  neglected  as  far  as  band  spreading  is  concerned. (  The  doppler  shifts 
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dc  however  effect  fading )  The  maximum  doppler  shift  Af^  is  computed  from 


<26) 


If  2°  and  f~  3GHz 


Afj  =  -  f  =*  (.  54)u  Hz  where  u  is  the  irregularity  velocity  in 


m 


/  sec. 


(d)  The  fading  rates  are  high  both  due  to  the  long  distance  and  high 

frequency.  The  mean  fading  time  is  of  the  order  of  0.  1  sec.  or  less. 

Receiver-  The  receiver  is  of  course  the  key  to  the  feasibility  of  this 
troposcatter  system.  Let  us  assume  that  we  are  going  to  operate  in  a  non 
coherent  FSK  mode  with  a  binary  bit  transmission  rate  of  1000  bits/sec. 

Let  us  assume  further  that  at  present  we  are  restricted  to  frequency  and 
time  diversity.  We  will  propose  the  following  receiver  system  so  as  to 
obtain  a  sizing  of  the  transmitter  power  (alternate  and  more  detailed  specifi¬ 
cations  of  receiver  systems  will  be  examined  in  a  future  report).  Each  bit 
is  constructed  of  three  frequency  chips  of  duration  l/3  millisec.  This 
specifies  6  possible  transmission  frequencies.  The  received  signals  at  the 
three  frequencies  that  constitute  a  bit  are  assumed  to  be  statistically  independent 
This  requires  that  the  three  frequencies  be  separated  by  the  propagation  channel 
bandwidth  mentioned  above.  Let  us  assume  that  we  have  an  optimum  (frequency) 
diversity  combining  system  based  on  Rayleight  envelope  statistics.  A  typical 
receiver  channel  is  shown  in  Fig.  16*  Zero  or  one  decisions  are  based  on  a 
comparison  of  the  values  of  the  y.  values.  [  12]  .  The  probability  of  error- 

versus  the  energy  per  bit  to  noise  ratio  is  shown  in  Fig.  17  with  the  number 

-  3 

of  chips  per  bit  as  a  parameter.  Then  for  a  bit  error  probability  of  10  and 
for  3  chips  per  bit  the  ratio  E/N  =  ( —  17.4  dB). 

Pm  e 

The  signal  to  noise  ratio  is  given  ^  w^ere  T  is  bit  duration 

(1  millisec.  )  Thus  the  signal  to  noise  ratio  in  dB  is 
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|l>-  6-9591 1 

TRANSMITTED  SIGNAL  s^tlcosoi,! 

RECEIVED  SIGNAL  *(t)  *  o(t)$0(«-r)  cos[w,(t -r)  +  $(t)]  +  n(t> 


Fig.  16.  Typical  receiver  channel 


PROBABILITY  OF  ERROR 


Fig.  17.  Binary  system  with  optimum  diversity  combining. 
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17.  4  +  30  =  47.  4  dB 


P 

m 

~w 

The  minimum  average  received  signal  level  for  this  receiver  system  must  be 

=  N  +  E/N  +  l/T  (terms  in  dB)  (27) 

P  [dB]  =  47.  4  -  203  =  -155.  6  dBw 
m 

This  is  the  value  which  will  be  used  in  the  transmitter  power  calculation. 

It  is  assumed  that  individual  receiver  bandwidths  need  only  be  of  the  order 

B ~  n  -L~  3000Hz. 

P  I 

We  have  in  a  sense  used  only  frequency  diversity.  One  disadvantage 
of  this  system  is  the  following.  Assume  the  propagation  channel  bandwidth 
is  1  MHz.  Then  the  total  RF  bandwidth  must  be  greater  than  Z  MHz  (for 
3  frequencies).  This  could  be  in  conflict  with  frequency  allocation  space  in 
the  GHz  range.  One  could  easily  envision  an  equivalent  time  diversity  system 
with  the  same  signal  to  noise  ratio  but  a  data  rate  reduced  by  a  factor  of 
3.  (the  total  RF  bandwidth  here  would  be  only  of  the  order  of  10  KHz).  To 
keep  the  logistics  simple,  spaced  antenna,  angle  of  arrival,  polarization 
and  multipath  diversity  have  not  been  considered.  A  combination  of  these 
methods  could  of  course  provide  economy  of  transmitter  power  at  the 
expense  of  complicating  the  operational  details 

I-  Sample  Transmitter  Sizing  Calculations 

We  will  assume  that  the  scattering  angle  is  6.  Z°  and  the  antenna 
beamwidth  ~  1.  5°.  Then 
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f  =  3  GHz 


Minimum  Received  Signal 


N  =  kTop 
Eb  ^NT 

Antenra  Parameters 

-203  dBw 

+47.  4  dB 

Gain 

-84.  6 

(d  ~4  +  , 
meters) 

Antenna  LTL 

Effi  ciency 

+  6.  0 

Antenna- to -Medium 
Coupling 

Free  Space  Power  Loss 

+  9.  0 

L  t 
pf 

Propagation  Channel  Loss 

+  158.  2 

Gaseous  and  Rain  ^G,  R 
Attenuation 

+3.  2 

Cross  Section  a 

Distance  Dependence  - 

-13.  6 

d2 

4rr lx '  |  2  Ifl-X'  |2 

+  1 15. 2 

Required  Transmitter  Powers  +37.  8  dBw 

or  6  KW 


=  4  GHz 

f  = 

5  GHz 

-203 

-203 

-203 

+47.  4 

+47.  4 

+47.  4 

-84.  6 

-85.  7 

-92.  6 

( d  ~  3  . 

meters) 

(d  ~2.  6, 
meters) 

(d  ~4.  0  . 
meters) 

+  6.  0 

+  6.  0 

+  6.  0 

+  9.  0 

+  9.  0 

+  14.  0 

160.  7 

162.  6 

162.  6 

+4.  2 

+5.  7 

+  5.  7 

-  14.  0 

-14.  35 

-  14.  25 

+  1 15.  2 

+  115.  2 

+  1  15.  2 

+40.  9 

+42.  85 

+41.0 

12.  3  KW 

19.  3  KW 

12.  0  KW 

It  is  seen  (with  the  exception  of  the  2.  6m  antenna  at  5  GHz)  all  the 
computations  give  values  in  the  neighborhood  of  10  KW.  We  feel  that  these 
values  are  rather  conservative  values  for  the  following  reasons. 

(a)  The  cross  section  <j  was  derived  on  the  basis  that  only  the  tropopause 
layer  contributed  to  the  scattering.  Obviously  both  the  'weaker'  stratospheric 
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and  upper  troposphere  layers  will  contribute  significantly. 

(b)  There  may  be  no  need  for  such  a  high  data  rate  (1000  bits/sec. 
normal).  Hence  we  could  either  reduce  the  bit  rate  (to  as  low  100  bits/sec. 
and  10  bits/sec)  and  pick  up  as  much  as  20  dB  or  we  could  increase  the 
probability  of  error  from  10“^  to  10-2  and  save  4.  7  dB  (or  we  could  do  both). 
Certainly  we  could  adopt  other  diversity  schemes  provided  that  they  are 
consistent  with  the  operational  requirements  desired. 
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Ill .  R  E  C  O  MM  END  A  T  IONS 


We  have  shown  that  a  mobile  troposcatter  system  with  a 
moderately  long  link  range  (640  km)  and  an  information  bit  rate  of  the  order 
of  1000  bits /sec  is  feasible.  This  is  a  rather  uncommon  type  of  tropo- 
s catter  system  . 


We  advocate  further  research  being  done  in  tropospheric  scattering 
from  the  tropopause  and  stratospheric  regions  (particularly  in  the  2  to  8  GHz 
region).  It  is  felt  that  the  results  from  scatter  experiments  in  this  area 
will  provide  the  fundamental  basis  for  a  class  of  operationally  flexible 
troposcatter  systems. 
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APPENDIX  A 

Essential  Equations  of  Scattering  Theory 


The  following  few  paragraphs  outline  the  derivation  of  the  equation  for 

the  power  scattered  by  the  refractive  index  inhomogeneities  in  the  common 

volume  VD  .  formed  by  the  intersection  of  two  narrow  antenna  beams.  If 
K 

^  <co  c  the  wave  equation  can  be  written  as 
o  t 

k  2 

VxVXE(r)-k2E(r)=k2  -  1)  E(r  )  —  P(r  )  (a-1) 

-*  -»  O  — *  — »  Of  f  -»  -♦ 

o  ^  o 

A  small  perturbation  technique  (the  Born  approximation)  is  used  to  solve 
this  equation.  The  electric  field  is  broken  into  separate  incident  and 
scattered  components  and  the  1st  order  equation  becomes 


V  X  v  x  E  SC(r  )  -  k  2ESC(r  ) 

— »  — »  O  — >  * 


o  if  (r  i  VRI) 
k2 

PlnCld-  (r,  r')  if  (r€  Vr,) 


(a-2) 


The  standard  method  of  solving  the  above  equation  is  to  use  the  well  known 
Green's  identity,  which  in  turn  first  requires  finding  the  simpler  Green's 
dyadic  solution  .  Now  to  solve  the  equations  for  the  Green's  dyadic,  one 
uses  Fourier  transform  calculus  along  with  the  usual  complex  function 
theoretic  arguments.  If  only  the  far  field  terms  are  included  the  scattered 
field  turns  out  to  be 


The  equations  for  the  Green's  dyadic  (and  its  compatibility  equation)  are 
given  in  Cartisian  tensor  notation. 


dX.dX.  dX  bX 
1  J  P  P 


-  k 


F  ■  =6.  u  (r  f,  r  ,  t,  t ') 

A  lm  lm  ov  g  •  s' 


ko  dX.  rim'  '6imaX.Uo(rg’rg',t,t'1 
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kV*0* 


T-,  sc,  ,,  o 

E.  r  ,  t)  =  - 7z- 

l  g  477 


I 

V 


R1 


6. 


lm 


r  -  r 


(r.  -r.')(r  -r  ') 
l  i  m  m 


r  -  r' 


e(r  ’) 

— s —  1 


jk  |r  -  r  '  I 

rpinc  ,  J  o  I -  1  ,  .. 

E  (r  ')  e  d  V 

m  g 


(a-3) 


R 


whe  re 


|r  -  r'  |  =  v  (rj  -  rj  'Hrj  -  rj ')  • 


s  c 

In  equation  a-3,  r  refers  to  the  field  point  where  E  is  measured  and  r’ 

in  c 

refers  to  the  scattering  source  point.  It  is  clear  that  is  an  inconvenient 

variable  and  should  be  related  to  the  field  near  the  transmitter  or  better  yet 
to  the  total  power  radiated.  The  geometry  to  be  used  and  the  symbols  are 
shown  in  Fig.  18.  The  transmitter  (A)  is  the  origin  of  the  coordinate  system. 
Using  the  notation  of  Fig.  18,  the  electric  field  at  the  receiver  can  be 
expressed  as 


E(d,t)  = 


k02e'JWt  6e(x')  aD(kr)sin  X  |E(0)| 


4n 


I 


x ' 


R 


jf  (k  -  kr)  .  x'  +  kr  .  d  } 

X  e  dV 


R 


(a- 4) 


where  ,  is  the  common  volume  and  |E(0)  |  is  the  magnitude  of  the  field 
at  the  transmitter.  Since  the  beams  are  narrow  (~  1°)  one  can  make  the 
approximations 


a)  k^,  =  k^,  and  kr  =  k for  all  x1 

b)  Gr*kx'^  =  GT*kX'^  GR*kr'  =  GR*kR^ 


(a-5) 


The  field  at  the  receiver  feed  E^(d,t)  is  then  related  to  the  field  at  the 
transmitter  feed  E^O)  by 
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Er  (d.t) 


aD  sin  X  (ET(0)  |kQZ/G^G^  \  -j(wt-kR.  d) 
- ; - -  e 

(417)  |X  '  |  |d-X'  |  / 4F 


6c(x')  j(kx,-kR).x' 

X  - * -  e  "*  d  x 1 


(a-6) 


5'c 

The  usual  practice  is  then  to  form  E^.  and  normalize  by  (E  .  E 

which  is  proportional  to  the  power  that  would  be  received  on  a  line-of- sight 
transmission  between  the  antennas.  The  free  space  power  is  proportional  to 


(E.e\  e  |ET(0)|2GTGR(A3)2 


T  R  '477  d 

Now  =  2  N  x  10  ^  and  hence  the  final  power  equation  reduces  to 


(a-7) 


P(^>*  _  5r*  Sr  _  sin2  X  K4d2  10" 12 _ 

Pfs  (E.  E')fs  1677  sin4  •&  |X'  |2  |d- X'  |2 

+  j  K.  (x'  -  x") 


J'  J  N(x')  X 

V  V f 


(a-8) 


N  (xn)e 


dx  'dx 1 


where  K  =  k  -  k  ~  and  K  =  2k  sin 

—  — x  -»R  o  2 


It  must  be  recognized  that  equation  (a-8)  is  appropriate  for  a  particular 
spatial  distribution  of  the  refractive  index  and  accordingly  both  N(x;t)  and 
P(d)  are  random  variables.  As  yet  the  structure  of  turbulence  in  the  common 
volume  (V^f)  has  not  been  specified.  Let  us  consider  the  term 

-jK  .  (x'-  x  "  ) 

J*  f  N(x»)N(x")e  dx»  dx'  . 

V’  V" 

If  N  (x1)  is  a  constant  this  expression  is  zero  for  K  ^  0  (which  one  expects 
physically).  A  cumulative  distribution  function  F^(x'  ,  xM)  is  assumed  and 
one  then  determines  the  expectation 
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-j  K  .  (x 1  -  x") 

f  J"  E  [  N(x' )N(x")]  e  dx"dx' 


+j  K.  (x'  -  x") 
e 


X  l  cov  N(x')N(x")  +  EN(x')EN(x")  (  dx"  dx' 

j  — »  — ♦  — ♦  — *  (  -4  — * 


(a-9) 


The  second  term  on  the  right  side  of  equation  (a-9)  becomes  zero  under  the 
integration  process.  If  the  field  is  locally  spatially  homogeneous  (in  the 
wide  sense)  one  can  write  cov[  N(?c ')N(xJ')]  =  Rj^(x x")  where  R^(’ ,*)  is  the 
correlation  coefficient. 


cov  (N(x')N(x")]  =  Rn(x'-x",  0)  =  Rn(x'-x") 


(a-  10) 


Now  we  perform  a  change  of  variables  (x x  ")  -♦(£  ,  y )  where  £  =  x 1  -  x  "  and 
x 1  +  x "  “* 

y  -  — —  .  The  Jacobian  for  the  change  of  coordinates  turns  out  to  be 

equal  to  1  hence; 

+jK.(x'-x") 

P  J  E  [  N(x')  N(x")]  e  dx'  dx"  =  j  [  P 


V  „  V  , 
x  "  x  1 


V  „  V, 
x "  £ 


+jK.  £ 

*  rn  e  d f]  dx  • 

The  normalized  autocorrelation  coefficient  pj^(«  )  is  defined  by 


Rn(£)  *  E|N2]  pN(|)  *  (N2)  pN(|) 

Y 

where  (N  )  can  have  slowly  varying  spatial  dependence.  A  spectrum  cf^.(K) 
of  the  autocorrelation  function  is  defined  by 

i  ,  jg-f 

~dS  =  <a 

8  ff  “* 

where  it  is  assumed  that  p^  decreases  rapidly  with  |£  |.  Therefore 
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-j  K  .  (x '  -  x  ")  , 

f  f  E[N(x')N(x")]  e  "*  dx "  dx '  =  (2  77 )  d>(K)  f 

' '  1  “ ♦  — »  1  — ♦  —4  — ♦  iJ 

V 


X  (N  ) dy  . 

Therefore  equation  a-8  can  be  written  as 

J  <N2)  <■* 


fs 


|X'|2|d-X’|2 


(2  sin  ) 


(a- 12) 


(a-  13) 


If  the  turbulence  has  the  same  intensity  throughout  the  common  volume  then 


$N(K)(hT)  V  (a- 14) 


P(d) 

~  o  2 

10"12d2 

4  2 

K*  sin  X 

p  —  877 

fs 

|X-|2  |d-X'|2 

(2  sin  -&-)4 

The  equation  for  bistatic  scattering  if  path  and  system  losses  are  ignored  is 


P(d) 


PtGtGr>  a 


(4ir)3  |X»  |2  |d-X'|2 


(a-  15) 


Normalized  by  the  free  space  power  we  get 

,2 


P(d)  _ 
Pfs  " 


ad 


IX')2  |d^-X'  |2  (4ff  ) 


(a-  16) 


If  equation  a- 16  and  a- 14  are  combined,  the  bistatic  cross  section  turns  out 
to  be 

K  4  •  ^ .  (a-  17) 


a  =  4(277-  )3  (10"12) 


(2  sin  & 


sin2  X  (N2)V$n(K) 


2  77 

Since  k  = 

o  X 


K 


2  .inf 


the  above  equation  can  be  written  as 


a  =  (4(217)7,  10'  12)  (N2)  V$n(K) 

X 
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or 


f  (2.  10'6)2  (N2)}  8tt  3  k4  sin2  X  V*n(K) 


This  result  came  from  assuming  that  the  turbulence  was  uniform  over  the 
volume  V  and  is  the  most  common  formula  for  the  scattering  cross  section. 
We  will  use  a  slightly  different  version  which  is  consistent  with  that  of 
Tartarski.  Instead  of  N  units  we  will  use  the  refractive  index  77  and  re¬ 
write  equation  a- 8  as 


and 


P(d) 


2  4  2 

sin  X  K  d^ 


3fs  (4)(l6ff)  sin4|-  |X'|2  |d -X1  |‘ 


f  f 

V  ,  V 


x  ’  x  " 


+j  K  .  (x1-  x") 

X  E  [  77  (x  ’)  77  (x ")  ]  e  dx" dxf 


(a-  18) 


and  let 


J*  f  E[p(x')v(x")]  e 


j  K  .  (x 1  -  x  ")  — - 

dx "  dx  1  =  J  77  (if)  J* 


V  ,  V  , 
X  1  x' 


y 


K 


jK.(f)  «  3 

X  PN(£  ^  e  d  £  d^  =  8V  f  <J>  (K,  y)dy . 

V 

y 


(a-  19) 


Then 


P(d) 


2  .  2  v  VA  .2 

n  sin  X  K  d 


fs 


8  sin4f  |X’|2  |a-X'|2  '  n 


/  (a-20) 


a  = 


■E-  sin  X  K4  r  Q  (K,v)dy 

2  Sin4£  J.  %  -  %  Z 


0  2  .2  ,  4  r* 

8n  sin  X  k  j 

V 


x 


(a-21) 
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In  terms  of  dB  the  appropriate  expression  to  be  used  is  then 


P(d)  =  P  [dB]  +  10  log  - - j  +  o'  [dB]  . 

4n  li"  X'f  |  X'p 


(a-22) 
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APPENDIX  B 

Approximate  Geometry  of  the  Common  Volume 


A  few  useful  but  approximate  formulas  for  the  common  volume  (for 
narrow  beam  antennas)  are  given  below.  The  beams  are  assumed  to  be 
rectangular  in  the  plane  transverse  to  the  beam  axis.  The  area  of  face 
ABCD  is  given  by 


A(ABCD) 


(dr  0V,  r)(dt  6V,t) 

sin  p 


The  width  of  the  common  volume  is  determined  primarily  by  the  narrowest 
beam.  Define 


W  -  d0H  “  mln  (dr  0H>  r’dt  ®H,t*' 


The  volume  of  the  common  volume  is  approximately 
V  -  de. 


^dr  0V,  r'1  (dt  0V,  t* 


H  sin  p 

If  the  geometry  shows  perfect  symmetry  d  =  D/2;  0  = 

=  0, 


V»  R 


V 


3  A3 

D  el 

8  sin  & 


(b-  1 ) 


0V,  R  0H,  R  0V,T 


(b-2 ) 


Another  useful  formula  is  the  area  of  a  horizontal  layer  which  is 
limited  by  the  common  volume  geometry. 

(AZ)D  0. 

±Ll  ±L  n  - 

n  *  ? 

tan- 


A  s.  o 

L  ^  *  2  0H 

2 


'H 


(tan  f) 


(b-3) 
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ll8-  6-9594| 


r 
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For  any  two  ray  paths,  the  maximum  path  difference  is  given  by 

Al 

For  the  case  of  beam  geometry  symmetry 


.  dr9v»r>  +  (dr6v^  -  dtflv.t) 

sin  ft  tan  sin  ft  tan  ft 


Al  =“  D  0  tan 


(b-4) 
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tropopause  layer  is  to  be  included  in  the  common  volume  of  the  beams.  This 
effectively  limits  the  range  in  most  cases  to  about  450  miles.  It  is  proposed  to 
use  frequencies  in  the  interval  3  to  5 GHz.  This  enables  one  to  specify  a  com¬ 
pletely  mobile  troposcatter  system  since  the  antenna  paraboloids  have  diameters 
of  the  order  of  10  feet  and  the  power  requirement  is  a  modest  10  kW. 
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